Dynamical role of anyonic excitation statistics in rapidly rotating Bose gases 
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We show that for rotating harmonically trapped Bose gases in a fractional quantum Hall state, the 
anyonic excitation statistics in the rotating gas can effectively play a dynamical role. For particular 
values of the two-dimensional coupling constant g — —2irh 2 (2k — l)/m, where A: is a positive integer, 
the system becomes a noninteracting gas of anyons, with exactly obtainable solutions satisfying 
BogomoPnyi self-dual order parameter equations. Attractive Bose gases under rapid rotation thus 
can be stabilized in the thermodynamic limit due to the anyonic statistics of their quasiparticle 
excitations. 
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Introduction. -Classical and quantum properties of 
atomic Bose gases when they are set under rotation 
were intensely studied in the last couple of years both 
from the experimental El 0, 01 and the theoretical sides 

HHSHHUHEJ ° f P articular interest is the 
behavior of rapidly rotating two-dimensional (2D) gases 

in a fractional quantum Hall state, corresponding to an 
electrically neutral, bosonic analog of ultrapure 2D elec- 
tron gases in very strong magnetic fields, which exhibits 
strongly correlated physics 0, 0| ■ It is well established 
that fractional quantum Hall states possess an yon ic ex- 
citations above an incompressible ground state [13, car- 
rying certain fractions v of the "elementary" charge , 
where v < 1 is the filling factor of the lowest Landau 
level. On the low-energy level of an effective field the- 
ory, these systems are described by their coupling to a 
fictitious statistical gauge field A 1 " 1 . This is known as 
a Chern-Simons effective field theory description of the 
fractional quantum Hall effect [lj, |l5| . 

In the thermodynamic limit, a system of untrapped 
bosons with attractive interaction is unstable against col- 
lapse and can stabl y ex ist only for a finite trapped num- 
ber of particles \jj&U^. The purpose of the present con- 
tribution is to point out that a 2D harmonically trapped 
rotating Bose gas with attractive interactions can be sta- 
bilized by the statistical gauge field associated with the 
fractional charge of its quasiparticle excitations. The 
anyonic nature of the excitations can therefore play a 
dynamical role, in that it can compensate the nega- 
tive coupling constant associated to the interaction of 
atomic bosons. This dynamical role is due to the fact 
that the statistical "magnetic" field is proportional to 
the density of the system, resulting in particle-statistical 
flux composites. For a special value of the statistics 
parameter = i//[2w(l — v)] or, correspondingly, for 
given 0, at a particular value of the coupling constant 
g = — 2nh {v^ 1 — l)/m, the interaction can effectively 
be eliminated entirely and make the system behave as an 
interaction- free gas of anyons. In physical terms, the sta- 
bilization of the attractively interacting gas in the frac- 
tional quantum Hall state against the collapse to a sin- 



gular distribution may be understood by the existence of 
a 2D analog of Fermi pressure, "anyonic" pressure. 

The order parameter distribution for the free anyon 
gas can be found analytically provided the Bogomol'nyi 
self-duality identities are satisfied yjj, which then leads 
to topological vortex solutions for the Ginzburg-Landau 
order parameter 

mm Hummus 

Rotating atomic 

Bose gases offer a unique opportunity to actually prepare 
and detect them rather directly and with great accuracy, 
and to study their breather soliton dynamics [2(1 l2lj |. 
Starting from the statistical mechanics of particles in the 
plane, interacting by contact potentials, other investiga- 
tions on obtaining a free anyon gas can be found in |25j ; 
however, only repulsive interactions were treated, with- 
out externally imposed magnetic or rotation fields. 

Abelian statistical gauge field.— We assume that we have 
a 2D gas in the fractional quantum Hall regime which 
admits a Ginzburg-Landau type description in terms of 
a Chern-Simons theory [l^. A statistical (2+l)D gauge 
three-potential A^ may be implemented by showing the 
physical equivalence of the two Hamiltonian theories with 
A>* = and A» 7^ on the quantum level Q. The 
primary effect of the statistical vector potential is that it 
allows for an additional Chern-Simons term in the low- 
energy action (Ti = 1): 



Se = 



[ d 2 xdt ( m*(d t + iAo)^ - — \DV\ 

J V 2m 



-V(x)\*\ 2 Um 2 ) + f e^VU^ 7 ) (1) 

Here, D is the gauge covariant spatial derivative defined 
in Eq. J5J, the Abelian field strength of the statistical 
gauge field reads T a p = d a Ap — dpA a , and W(|^| 2 ) is 
the interaction energy density. We use relativistic no- 
tation, i.e., a,/3,7,/i = 0,1,2; the spacetime metric is 
diag (1, — 1, — 1) and repeated indices are summed over. 
All higher derivative terms in the statistical gauge po- 
tentials (e.g., Maxwell terms oc FafiF ) are left out in 
this low-energy, low- momentum expression |2^ . 

Classically, the momenta of the particles making up the 
system are given by the decomposition p = mv + nA+A, 



2 



(the spatial part of A^ is denoted A), so that the gauge 
covariant derivative in the quantum domain is given by 

(2) 



D = V ikA - iA. 



Here, A is the external U(l) vector potential of rotation, 
yielding the Coriolis force, with V x A = 20, twice the 
applied rotation field f2 perpendicular to the 2D plane, 
and k — m is the coupling constant for a rotating Bose 
gas. The total scalar potential 



V 



Vtrap - -mtt 2 r 2 



(3) 



consists of (harmonic) trapping potential and centrifu- 
gal potential. Due to the Chern-Simons term, the field 
strength of the statistical gauge field is given by the cur- 
rent as follows: 
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(4) 



The homogeneous Maxwell equation for the statistical 
field strength, d^T a p + dpT^a + daTp^ — 0, is then au- 
tomatically equivalent to current conservation, 9 M J M = 0. 

The statistical gauge field strength is according to (0J 
dual to the current. This means, in particular, that the 
statistical magnetic field is proportional to the density: 



-QB = |*| 2 =p. 



(5) 



That is, the order parameter modulus is inextricably 
linked to the statistical flux, and particle-flux compos- 
ites are formed. Taking 8 > 0, by fixing Q. > 1, the sta- 
tistical magnetic field cancels part of the applied "mag- 
netic", i.e., rotation field in the canonical momentum p, 
and the vector potential effectively acting on the par- 
ticles is reduced. Defining the Landau level filling fac- 
tor of the original bosonic particles to be v = Trpo/m^t 
[3. ITU, with po = | 2 a homogeneous background den- 
sity, we have O = ^/[27r(l — v)]. In particular, for the 
v = 1/2 anyons discussed in yjj, O = 1/2%. The effec- 
tive magnetic flux of particle-statistical flux composites 
$ = § dx- [A+A/m] — ^<&o, then is reduced compared to 
the bare rotational flux quantum $o = § dx ■ A = 2ir /m 
obtained for vanishing A |27| . 

The Ginzburg-Landau energy functional of the gas in 
the rotating frame, corresponding to is composed as 
usual of kinetic, scalar potential and interaction energy, 



H = 



J d2x (i 1 ^ |2 + ^wi*! 2 +^(i*i 2 ) 



(6) 



In the rapid rotation limit f2 ~ oj±, Vtrap = ^muj 2 _r 2 is 
(nearly) cancelled by the centrifugal potential, V ~ 0, 
where ui± is the trapping frequency perpendicular to the 
axis of rotation. 

Using the Bogomol'nyi decomposition for the ki- 
netic energy term in the integral 



and relation JSJ, we can rewrite the Hamiltonian in two 
different ways corresponding to the ± sign in , to read 



H= d 2 x — \(D 1 ±iD 2 )^\ 2 + V(x)\V\ 
2m 

2 V mb 



(8) 



We have neglected the term involving the curl of the mat- 
ter current ±V x J, as it amounts, after integration, 
to a surface contribution vanishing for sufficiently well- 
behaved current fields and/or well outside the boundary 
of the rotating gas. We infer from the last term in the 
second line of © that due to the proportionality of den- 
sity and statistical magnetic field expressed by Eq. JSJ, 
and employing the identity 0, part of the kinetic en- 
ergy \D^\ 2 /2m can effectively act as interaction energy, 
The potential in the Ginzburg-Landau energy func- 
tional may generally be expanded 



^(l*| 2 ) = fl*| 4 



Z|$|6. 

6 1 1 



(9) 



where g is the value of the two-body coupling constant in 
the fractional quantum Hall state; for sixth order stabil- 
ity of the system, 7 > is required. Taking into account 
the coefficient of the quartic term in ® resulting from 
the above expansion and Eqs. (QJ and J7J), the effective 
interaction coupling of the system is now defined as 



ffoff = 9T 



1 

mQ 



(10) 



{D 1 ±iD 2 )^\ 2 ±V x J ±Bp±2mflp, (7) 



If we choose the relation of coupling constant and statis- 
tics parameter to be g = ±l/mO, we see that we have 
effectively eliminated the quartic interaction coupling 
term. By "interaction-free" we mean that we have can- 
celled the interaction g by a proper choice of the statis- 
tical magnetic field B, while the covariant derivatives D\ 
and Z?2 do of course still implicitly contain the statistical 
gauge fields A\ and A2, and therefore, due to (QJ, the 
density distribution. 

We conclude from relation 11U|) that a large origi- 
nal coupling g can be converted to a g c ff approaching 
zero. Due to the statistical interaction, a negative g 
system, unstable towards collapse in the thermodynamic 
limit 16], may be stabilized by the statistical interac- 
tion, in the sense that there are nonsingular distribu- 
tions "i>(x,t) (see below), with finite energy. The value 
9 = 9c = —l/m@ defines a critical negative coupling 
strength, below which the gas cannot be stabilized for a 
given O. 

In the limit that the coefficient of \^\ 4 vanishes, g c s — 
in (fin)) , the problem of determining the ground state 
of zero energy becomes exactly solvable (neglecting the 
small I "J 1 6 term in the Ginzburg-Landau expansion). The 
energy per particle then assumes its lower limit Hq /N = 
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±fi, provided the self-duality (Bogomornyi) constraints 

m 

{Di±iD 2 )$ = Q (11) 

are satisfied. To protect the self-dual states of effectively 
zero coupling, there is an energy barrier to neighboring 
states in ^ space, which do not have the property that 
9cS — 0. The magnitude of this energy barrier depends 
on the (positive) contribution of the \^\ 6 term in the 
Ginzburg-Landau type expansion of W(|^| 2 ). 

Jackiw and Pi [lj| have shown that solving the self- 
duality equations Ijlll) above, in the case of zero "exter- 
nal" field, i.e., in the problem with just the statistical 
gauge field present and trapping and rotation turned off, 

V = and Q = 0, is equivalent to solving the Liouville 
equation 

Amp = ±2/9/6 (V(x) = fi = 0), (12) 

whose complete set of solutions is known [28|. To ob- 
tain regular and non-negative solutions of the Liouville 
equation, the sign of the right-hand side must be cho- 
sen opposite to that of 0. Thus the lower, minus sign 
is appropriate for 8 > 0. The general solution of the 
Liouville Eq. I|12|) in the plane of complex z = x + iy 
reads p(z) = 46|/'(z)| 2 /[l + |/(z)| 2 ] 2 , where f(z) is an 
arbitrary holomorphic function. Radially symmetric vor- 
tex solutions, for which 'J = .y//jexp[m</>], take the form, 
choosing f(z) = (z a /z) n , 

46n 2 [/r \» (r \ n "T 2 

where n is an integer, and ro an arbitrary length scale 
reflecting the scale (dilation) invariance of (|12fl : for r — > 
0, p(r) oc r 2 ( n ^ 1 ) ; and for r — > 00, p(r) cx r - 2 (™+ 1 ). 

The solutions of (|1 1|) with a constant external mag- 
netic/rotation field (corresponding to our being in the 
rotating frame), and a linearly increasing electric field 
in the plane (corresponding to our linear trapping and 
centrifugal forces) can be obtained from the vortex soli- 
tons of Jackiw and Pi: Adding these additional external 
fields the problem (still) is quadratic H3, EL H3| The 
scalar potential in the Hamiltonian JSJl, i.e. the effective 
harmonic trapping field, is very small close to criticality, 

V = \muj\r 2 — imfJ 2 r 2 ~ 0, and can be neglected. The 
classical problem of finding the solution of (|llfl in an ex- 
ternal rotation field then corresponds to the problem of 
finding the (semiclassical) Landau levels of anyons [2Cj . 

To generate time-dependent vortex soliton solutions, 
one constructs a coordinate transformation, whose in- 
verse effectively removes the external field and thus leads 

us back to csj 23 mm. The most important feature, 

apart from the cyclotron motion executed by the soliton, 
which results from this transformation, is that the size 
of the soliton breathes if the background 7^ 0: The 



scale factor in (|13f) becomes time dependent according 
to ro — > ro cos [fit]. Hence the soliton size oscillates with 
the applied rotation frequency. Furthermore, the energy 
of the soliton diverges for n = 1 and finite ro, and thus 
according to (|13fl only n > 2 configurations with p = 
at the vortex line can be generated; note that for n = 1 
the density (|13l) at the origin has a finite value. The 
statistical vector potential A decreases at large distances 
from the center of the vortex line like 1/r, as required for 
a topological vortex of a given quantized circulation. It 
should also be stressed that the feature that there exists 
a breather soliton solution in the presence of the exter- 
nal rotation field, persists when g g c and there is no 
self-duality fulfilled according to lllfl : however, in that 
case simple analytical solutions like the one displayed in 
(|13[) cannot be obtained, because the Ginzburg-Landau 
equations then remain essentially nonlinear. 

Rescaling the coordinate vector via x = x^Q/2p , the 
typical length scale of inhomogeneous solutions of the 
equations l(TT|l is set by £0 = \/Q/2po = y/l/2mpo\g\, 
the analog of the coherence length in the repulsive case. 
The topological, i.e., quantized circulation vortex solu- 
tions in a rapidly rotating gas with g < 0, resulting from 
the boosted solutions (|13[1 , then are cousins of their coun- 
terparts in the repulsive-interaction superfluid. In the 
latter case, the core size of the vortices depends on the 
applied rotation rate |9| (for an experimental verification 
see the second reference of Q). Here, by contrast, the 
scale ro in Eq. (| 1 31) is essentially a free parameter, because 
the total energy of the soliton is minimized for r = (2l| . 
Therefore, an external force, created by a blue-detuned 
laser for example, has to be applied to the gas to gener- 
ate a soliton with a finite value of ro . It should be noted 
that 11 ^ is strictly necessary to obtain any solution 
of nonzero V&O; i-e., the symmetry breaking expressed by 
"to 7^ is induced by the rotation of the gas: The solu- 
tions of i|12[). e.g., the vortex solution in Eq. (|13|l . are all 
asymptotically vanishing, = if SI = 0. 

Non-Abelian case. -The previous considerations can be 
generalized to the case of a noncommuting statistical 
gauge field (cf., e.g., [13,0]), which is of potential rel- 
evance in rapidly rotating spinor gases. An outline of 
such a non-Abelian generalization reads as follows. The 
non-Abelian field strength is 

T aS 3 = d a Ap - dpAa + [A a ,A@] (14) 

where the noncommuting gauge field reads A^ = A^T a , 
using the Lie algebra of the anti-Hermitian generators T a 
of the non-Abelian group, [T a ,T b ] = f a bcT c (summation 
over a, b, c is implied; f a bc are the structure constants of 
the Lie algebra). The Chern-Simons-Gauss law (JSJ reads 
in its non-Abelian form 

B=^T a (^T a ip), (15) 

where ip is the order parameter field, transforming ac- 
cording to some given representation of the gauge group. 
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The scalar and vector interactions in a spinor conden- 
sate may be parametrized to read (2^, [3(| 

Mnt = J^tyJ^i " fWV^IjS^ , (16) 

where the summation indices i,j,k,l cover the matrix 
index of the particular representation of tp chosen. We 
infer, using (|15fl , that the effective coefficient of the quar- 
tic spin-spin interaction term may in analogy to I|1U|) be 
defined, restricting ourselves to the lower sign, 

(c 2 )cfr = c 2 + — . (17) 

The statistical interaction therefore is able to change, as- 
suming co < 0, the spin-spin interaction from c 2 < 
(ferromagnetic) to c 2 > ("polar" 29] ). However, for 
a two-parameter interaction Hamiltonian like (|16|l . with 
Co 7^ 0, there is no self-dual Bogomol'nyi point in pa- 
rameter space. Only at Co = 0, such a point exists; 
then, the system, in general, can support vortices obey- 
ing non-Abelian fractional statistics. The matter den- 
sity components associated to these fractional vortices 
are solutions of the Toda equation, which generalizes the 
Liouville Eq. JEJ jiiE3|. 

Conclusion.-We have established the fact that a 
rapidly rotating Bose gas with attractive self-interaction 
can be made manifestly stable in a given fractional quan- 
tum Hall state, such that for certain values of the nega- 
tive coupling constant the Ginzburg-Landau matter wave 
field exhibits self-duality. The resulting free anyon gas 
interpolates between the extreme cases of bosons (zero 
statistical pressure at T = 0) and fermions (maximal sta- 
tistical pressure at T — 0) , and the original bosonic gas is 
protected against collapse because of this effective "any- 
onic" pressure. The preparation and experimental veri- 
fication of this particular fractional quantum Hall state 
should, in principle, be possible starting from a rotat- 
ing Bose gas with small positive interaction coupling and 
high angular momentum L = N(N - l)/2v = 0(N 2 ). 
To avoid the problem of stabilizing a bulk gas at these 
very large angular momenta, one can conceive of putting 
the system on an optical lattice |31| . Starting from the 
small positive interaction coupling Bose gas, one switches 
over a close-lying Feshbach resonance to the negative 
coupling strength side of the resonance. At the end 
point of the sweep a coupling is chosen which fulfills 
g = — 2ir(2k — l)/m, corresponding to a fractional quan- 
tum Hall state at filling v = l/2fc. 

I gratefully acknowledge discussions with Petr O. 
Fedichev and Nils Schopohl. 
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